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Abstract. The massive nonsymmetric gravitational theory is shown to possess 
, a linearization instability at purely GR field configurations, disallowing the use 

of the linear approximation in these situations. It is also shown that arbitrarily 
small antisymmetric sector Cauchy data leads to singular evolution unless an ad hoc 
condition is imposed on the initial data hypersurface. 

h=s : 

oo 

(N ■ PACS numbers: 04.50. +h, 04.25.Nx, 04.20.Ex, ll.10.Jj 

> 

o 

CO The Nonsymmetric Gravitational Theory (NGT) (Moffat 1979, 1984, 1991) grew out of 

a reinterpretation of Einstein's unified field theory (Einstein 1945; Einstein and Straus 
1946) as a theory of gravitation alone, thereby bypassing the problems associated with 
the interpretation of the additional structure as representing the electromagnetic field 
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(Infold 1950, Callaway 1953). There were indications of problems even with this 
I , reinterpretation (Kunstatter et al 1984; Kelly 1991, 1992) which were clarified by 

Damour, Deser and McCarthy (1992, 1993) who showed that the wave solutions of 
the weak field equations did not decrease at large distances from the source along the 
forward light cone (discussed in more detail in Clayton (1996a)). 

The theory was subsequently altered into what will be referred to as the massive 
Nonsymmetric Gravitational Theory (mNGT) (Legare and Moffat 1995; Moffat 1995a, 
b) by requiring that the linearized field equations reduce to those of a massive 
Kalb-Ramond (1974) field, guaranteeing that the linearized fields are well-behaved 
asymptotically far from the source (Clayton 1996a). The action for the Massive theory 
is (16ttG = c = 1) 



Sngt — 



J d 4 xE[~g AB R A s B - g AB W e[A [W] B] + l A A A 

+ l g ( AB )W A W B + 1 I m 2 S [AB] g[AB ] ], (1) 



where the first three terms are identical to the vacuum UFT action, the last term is 
the Bonnor (1954) term, and the remaining contribution is the new feature of mNGT 
designed to obtain the correct linearisation. (Note that the Lagrange multiplier l A is 
not a new feature in the NGT action; if one uses the connection as defined by Moffat 
(1991), it would not appear.) 
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The convention of this work is to define the inverse of the (in general nonsymmetric: 
gAB ^ gBA) fundamental tensor by g AB gBC = gCBg BA — Sq. Densities with respect 
to the determinant of the fundamental tensor g :— dct[gAB] will be indicated by 
boldface: T AB := ^/^gT AB , and the symmetric and antisymmetric part of a tensor 
will be indicated by T(ab) ■= (Tab+Tba)/2 and T[ A b] ■= {Tab — Tba)/2 respectively. 
A torsion-free covariant derivative is employed, characterised by the connection 
coefficients that define the parallel transport of the frame by V eA [e]s = r^ B ec (which 
is not required to be compatible with any tensor), and what are normally treated as the 
antisymmetric components of the connection coefficients are considered as a separate 
tensor A^ c , the trace of which is defined by := A B B . The action (1) has been 
written in a general linear frame, where the basis vectors are related to a coordinate 
basis via a vierbein ca = Ea^B^, generally leading to nontrivial structure constants 
[e A , e B ] = Cab C 'ec- 

The NGT Ricci tensor that appears in the action can be split up into two 
contributions: R^b = Rab + Rab': the nrst i s identified as the Ricci tensor (i.e., 
it reduces to the GR Ricci tensor in the limit of vanishing antisymmetric sector), and 
the second contains contributions from the antisymmetric tensor field Ag C 

Rab = ec[T c BA ] - ^b[T c ca ] - M r §c] + \^ac] 

i j^E "pC 1 j^E "pC 1 I t^E j^C 
1 BA l CE ~ 1 CB l EA^ 1 [AB] L EC-> 

R% = V ec [K] C AB + V e[A [A] B] + A^l D Ag C . 
In addition to the vanishing of the torsion tensor 

rpA . t^A r A r< A _ n 

J BC •— 1 BC ~ 1 CB ~ <^BC — U, 

the field equations of mNGT take on the form (Clayton 1996a): 
A A = 0, 

\ A = W AB) W B , 

Ve B {z] lAB] = & AB) W B , 

V ec [g] AB - K DB A A D g AD A B DC + ^ A 6 B V DE) W E = 0, 
Rab + V e[A [VK] B] - f WaW b - \m 2 M AB - 0, 
where the mass tensor appearing in the last of these is defined by 

MAS = g[AB] ~ gcAgBDg [CD] + \gBAg [CD] g[CD]- (4) 

By requiring that the linearisation of the antisymmetric sector take on the form of 
a Kalb-Ramond field, the field equations have picked up three additional constraints at 
the linearized level, which have since been shown to be absent in the nonperturbative 
theory (Clayton 1995, 1996b). The purpose of the present work is to examine the form 
of the exact field equations (3) near initially GR field configurations, demonstrating the 
fact that the exact evolution of the system is not well approximated by the linearized 
system. It will therefore have been shown that mNGT is not linearization stable about 
GR spacetimes, and so the linearization which guided the construction of the action 
does not describe the theory even in an approximate sense. Furthermore, it will be 
shown that unless an ad hoc constraint is imposed on the initial data hypersurface, 
unboundedly large velocities occur for arbitrarily small antisymmetric sector fields. 
This indicates a Cauchy instability (as defined in Hawking and Ellis (1973)), since the 
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field configuration at an infinitesimally small time later will not depend smoothly on 
the initial data near a GR configuration. 

The paper will consist of two main sections: the first will review the field equations 
of mNGT in 3 + 1 decomposed form, and the second will explore the above issues, 
making it clear how and why the naive linearization fails. The discussion of the surface 
decomposition of mNGT is fairly terse, and further details about the formalism are 
given in (Clayton 1995, 1996b). This construction is a straightforward extension of 
that given by Isenberg and Nestor (1980), although it must be admitted that the 
structure of NGT introduces a great deal of algebraic complication that necessitates 
the introduction of many auxiliary fields in order to present results in a reasonably 
compact form. In addition, in order to work consistently with a surface adapted basis 
in NGT, it was necessary to develop a 'moving frame' formalism for NGT (Clayton 
1996a). 

2. The Surface-Decomposed Field Equations 

What will be given here is a summary of the space and time decomposition that was 
used to consider the Cauchy problem in NGT (Clayton 1995, 1996b). Spacetime has 
been foliated into spacelike hypersurfaces (M,g) <~ {(E t ,g t )\t £ I C K}, where the 
metric of spacetime has been operationally identified with the symmetric components 
of the inverse of the fundamental tensor of NGT g( AB ~> . Configuration space has 
therefore been chosen to consist of g[ ±a l and g ab , where {a, b, . . .} e {1, 2, 3} indicate 
components of a tensor on £ expanded in a coordinate basis, and the index '_L' 
indicates a component normal to the hypersurface as defined by g( AB ). 

This surface compatible basis is designed so that the symmetric components of 
the inverse of the fundamental tensor take on the ADM form (Arnowitt et al 1959): 
g" 1 = e^®e^ — 7 < - ab -'e a <g>efc, where the basis may be written in terms of the coordinate 

basis (d t ,d a ) as: e± — j^dt — -%-<9 a , e a = d a , and the covector bases as: 9 1 - — Ndt, 
9 a = dx a + N a dt. The time vector is decomposed as t = Ne± + N, where (N, N a ) 
are the lapse function and shift vector respectively, and the notation ./V := N a e a has 
been employed. The non-vanishing structure constants are easily found to be: 

CL/ = e a [ln(A0L C ±a b = ^e a [N% (5) 

and the fundamental tensor and its inverse become: 

g = F9 ± ®9^-(a a ~f] a )e ± ®6 a -(a a +[3 a )e a ®9 ± -G ab 6 a ®d b , 

i h (M 

g = e ± <g> e ± + B a e± <g> e a - B a e a <g> e± - j ab e a <g> e b , 
where: 

F :=l + G (ab) B a B b , a a :=G [ab] B b , a := G (ab) B b , (6b) 

from which one finds that a a B a = 0. The spatial part of the fundamental tensor is 
given by 

( 7 ac - B a B c )G cb = G 6c ( 7 ca - B c B a ) = 5 a b . (6c) 

Note that invertibility of the fundamental tensor in (6 a) is required in order that the 
volume element in the action be nondegenerate, which in turn requires that G ab exist. 
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The spacetime covariant derivatives have been decomposed as usual (Isenberg and 
Nester 1980) into surface covariant derivatives (written as vi 3 -*) whose action on basis 
vectors in TS (the space of vectors tangent to £) is given by V a [e]& = rjj b e c , and 
contributions from derivatives off of £ which will be written in terms of the surface 
tensors: 



r :— r^j_, c a :— r oJ _, a a :— Tj^ a , a a :— Y\ j_, 

„,,a pa „,a po ; t-i_L 

W b 1 ±61 u b L bJ_> K ab ■— 1 a fc 



(7a) 



In GR, (7a) are related by the algebraic compatibility conditions: a a = 7 ab afc, c a = 0, 
r = 0, and u a b = j ac kb c , as well as the conditions that enforce the vanishing of torsion: 
a a = C± a , w a b = u a b + C±b a , r° fcc ] = 0, and fc[ afc ] = 0. The antisymmetric tensor A^ c 
is decomposed as: 

b a :=Ai±, Ja b :=AL v a b--=Kb, Kc--=Kc (7b) 

where by definition j a b is an antisymmetric surface tensor. The vector field Wa that 
appears in (1) is decomposed as Wa = (W, W a ), and the traceless part of X bc is defined 
by X^bc — Kc ~ &[bbc] ■ It is also useful to make the following definitions: 

K a b ■= k a b+ jab, (8 a) 

k% := \(ci ac K hc + 7 ca i^ c6 ) = l (ac) hc + l [ac hbc, 
j a b := l^Kbc - j ca K cb ) = 7 [acl fc 6c + 7 (ac) j6c 

the first of which is defined for notational convenience, and the remaining are 
symmetric and antisymmetric sector contributions from tensors that have had indices 
'raised' by j ab . 

In the Hamiltonian formalism of mNGT (Clayton 1995, 1996b) the canonical 
momenta are weakly equivalent to the densitized components of the fundamental 
tensor appearing in (6a), and the set of conjugate pairs of phase space coordinates are 



{{B a ,W a ) 1 h a \Kab)} 1 (9) 

where the field 

W a := -W a + 2b a , (10) 

naturally appears as as the momentum conjugate to B a . (In (Clayton 1995, 1996b) 
there is an additional pair of symmetric sector coordinates that does not exist in 
typical Cauchy analyses of GR which is removable by making imposing the second 
class constraints: p w ^/^g[B a , j ab ] , and u w ^ ab K a b strongly.) 

The fact that the spacetime covariant derivative has been defined to be torsion-free 
(3a) as in GR results in k^ a b\ — ^[ ab ] = ®> now combined with: 

a a ~c a + C^, w a b :=u a b + C ±b a , (11) 

and equation (36) gives: 



v\ « 0, b a « -X b ab . (12) 



The symmetric sector Lagrange multipliers that appear in (7) may be solved for 
in terms of the antisymmetric sector Lagrange multipliers and Cauchy data as: 

T^B a b a -\ 1 ^l3 a {W b -2b b ), 

(13a) 



c a w - Jab B b + ±/3 a W + \G [ab]1 {bc \w c - 2b c ), 
u% « B a 6 b + B c A" b + k a b , 

and using the solution for v a b given in (13c) one finds the surface compatibility 
condition that determines L bc in terms of Cauchy data: 

V (3) [7 ](a6) M 7 (ah) Cc + 2B ( V] + 7^ A^ c + 7^ A^. (136) 

(The Lagrange multipliers v a b must be replaced in (136) in order to solve for the 
surface connection coefficients, since v a b in (13c) has explicit dependence on r bc .) In 
the antisymmetric sector one can solve for: 

W « fvi 3 )[B]<\ v% « V b 3) [B] a + j a b - ±,5 b a W, (13c) 

leaving 

v (3)[ 7 ][ah] _ _2B[ a w b ] +7 [ab] c c + 7 (ad) A^ c + ~f (db) \ a cd + |4°7 (61e) (^e - 26 e ) (13d) 

to determine the Lagrange multipliers 6 a and A^, bc . 

The general solution of the full set of Lagrange multipliers is highly nontrivial, but 
will not be required here. Only the relation that determines b a in terms of Cauchy 
data (isolated by taking the trace of (13c?) and using the remainder of (13) in order 
to replace the Lagrange multipliers with Cauchy data), 

+ e b frW] + k b b B a - k a b B b + 7 [o6] ifec-B c - \l [ah] f3 b V^ [B] c , 

will be necessary in order to demonstrate the result. It is this rather surprising ability 
to isolate a relation that determines b a in terms of Cauchy data that allows the analysis 
in the next section to be performed in a fairly straightforward manner. 
It will be useful to define the two tensors: 

Of := 7^ - B a B b , Of := B a B b - \"f [ac] G [cd]1 {db \ (15a) 
the inverses of which will be denoted 0^ ab and 0^ ab respectively, given in general by: 



B a B b 



o-i-^- 



B a G b + G a B b B-B + j-G ^ {Wb) 
hab G-B + { 1 -B){G-B) Ualb 



Indices have been 'raised' or 'lowered' by ^ ab ^ and its inverse S ab (for example, 
G a := S ab G b , and ^ ab lS bc = 6%), and the notation B ■ G := B a G a has been 
employed. The 3-dimensional antisymmetric symbol e abc and the surface volume 
element v^S := V det S ab have been used to define: 

VS J (16) 



G[ab] = 2\ r Se abc G c \ G a := ^=e afcc G[ bc ] 



where e 12 3 = e 123 = +1. Clearly in the limit of vanishing antisymmetric sector, 0^ a \ 
in (156) is not in general well-behaved due to the presence of (j-G)^ 1 which is singular 
in that limit. 

Using these and the definition 

\V^,S a := e 6 [ 7 [a6] ] + k b b B a - k a b B b + j^ 3bc B c - ^(3 b V^ [B] c , (17) 

(14) is written as 

Ofb b « \OfW b + \OfW b + ±S Q , (18a) 
the solution of which is (in terms of the phase space coordinates) 

K « Wa + h°2a\O bC W c + \0^ b ~t>. (186) 

In terms of the configuration space variables, equation (14) becomes 

Ofh = \0?W b + \OfW h + ±S a , (18c) 

where W a must be solved for in terms of time derivatives of the configuration space 
variables (B a ,7 ah ) in (19a) below. The solution of (18c) may then be determined to 
be 

b a = \W a + \0~l b O\ c W c + \0~^\ (18d) 

resulting in 

W a = -W a + 2b a = 0- a \O b 2 c W c + 0- a \Z b . (18c) 

The evolution of the canonical momenta as determined from Hamilton's equations 
or, equivalcntly, from the compatibility conditions in (3e) that involve time derivatives 
of the fundamental tensor, is given by: 

B a « £f[B] a - iV(vi 3) [7] [afcl +7 [a6] C±6 ± + ll (ab) {W b - 26 b )), (19a) 
^(06) „ £g)[ 7 ](^) + N(j^(T + u) -2{-y (c(a '>u b l -~i^v b l), (196) 

^[ab] _ £(3)[ 7 ][a6] _ N (B a (7 b - BV - (T + u) - 2{^U b \ - J^V^), (19c) 

where £^ is the surface Lie derivative. Partial derivatives with respect to the time 
coordinate have been represented by an overdot (e.g. B a ), and are equivalent to the 
properly defined Lie derivative off of the surface (Iscnbcrg and Nester 1980; Clayton 
1995, 1996b). The evolution of the coordinates (from Hamilton's equations or (3f)) is 
determined to be: 

W a « £f[W] a - 2NZ [a±] + ±m 2 NM [a±h (20a) 
K ab « £f [K] ab - NZ ab + \m 2 NM ab , (206) 
where Zab as determined from (3/) are given by: 

Z(ab) ■= R\ ab) - v ( a Hb) + (r + u)k ab - a a a b - k cb u c a - k ac u c b 

+ b a b b - Jac v c b - j bc v c a + X ad X d bc ~Wa- 2b a ){W b - 2b b ), (21a) 
Z[a±] ■■= ^ a 3) [W] + \W{a a c a ) v£ 3) kf a 

- Jab<J b + (c b - a b )v\ + ub a - u\b b - u b c \ c ab , (216) 
Z W ■= -^[ai W U + 2V[ a 3) [6] fc] + V< 3 >[A]^ + (r + u)j ab + a c X c ab 

+ a a b b - a b b a - jcbU c a - j ac u c b - k ca v c b + k cb v c a . (21c) 
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In (21a) the surface Ricci tensor is denned from (2 a) to be 

R§! = ecKJ - W^la) + e a \Ftt) + n a rj c - r c da r d bc , (22) 

which is constructed from contractions of the intrinsic curvature of V^ 3 ^ on S. 

The remaining field equations in (3f) are the Gauss relation or Hamiltonian 
constraint (Clayton 1995, 1996b) 

U « -(«u - u%u\ + v%v» a + T^iZgij) + W{ P M^ + j ab M ab ) 

- 7 (Qb) (^Ta^Tbc ~ IW a W b + IW a b b ) + IV^g(W) 2 

- l lab] (-Vi 3 ) Wh + 2Vi 3 > [b] b + V< 3 > [X] c ab ) - 2uB% a 
+ Vi 3 ' [BV 6 + 7 M& 6 ] - c q (BV 6 + ^ a H b ) « 0, 

and the Codacci relations or momentum constraints 

W a « B h e a [W b ] - e b [B b W a ] + j bc e a [K bc } + V=ge a [u] - 2e b [k b a ] » 0. (236) 

Diffcomorphism invariance of the mNGT action guarantees that the Poisson brackets 
of the Hamiltonian and momentum constraints satisfy the same closing relations as GR 
(Tcitelboim 1973, 1980; Isenberg and Nester 1980), where the fact that ^( ab ) appears 
in the Poisson bracket of the Hamiltonian constraint is a consequence of the choice of 
using g( AS ) to define the unit normal to E, and is discussed in more detail in (Clayton 
1996b). 

One now has the exact evolution equations of mNGT, written in surface 
decomposed form that is suitable for considering the Cauchy problem in more detail. 
As a first-order system, one would first solve (13) for the Lagrange multipliers replacing 
them everywhere in Hamilton's equations (19) and (20) and the constraints (23). 
There would then be four undetermined Lagrange multipliers (N and iV a ) and four 
diffcomorphism constraints (23) on the Cauchy data (g). To consider the system 
as second-order, one would go further and solve (19) for the canonical coordinates 
(W a , K a b)i replacing them in the constraints (23)as well as (20), the collection of 
which would then become the second-order Euler-Lagrange equations of mNGT. 

3. Linearization Instability 

Considering the complexity of the field equations, it would be somewhat too optimistic 
to attempt to draw definite conclusions about global solutions of the mNGT field 
equations (see Moffat (1995c) for some results). Instead it will be sufficient to consider 
the instantaneous problem here, considering Cauchy data for which the antisymmetric 
sector will be chosen to be an arbitrarily small perturbation of the symmetric sector. 
Normally one would expect that that the perturbation of the initial data would lead 
to evolution of the system that may be considered (at least for small enough times) 
as a perturbation of GR evolution. All contributions from the antisymmetric sector 
to the symmetric sector should appear at second-order, and the evolution equations 
would not normally drive the perturbations to become large; indeed, this is precisely 
what the naive linearization of the system leads one to believe. 

Instead one finds in general that the exact field equations give contributions to the 
symmetric sector that appear at background order (i.e., not as a small correction, but 
of the same order as GR effects) , and the field equations that determine the evolution 



8 



of the perturbations result in arbitrarily large velocities. This not only invalidates 
the use of the linear approximation, but also explicitly shows a Cauchy instability in 
mNGT. 

To begin, we will show how the results of the linearizations given by Moffat 
(1995a, b) and Clayton (1996a) (about a fixed GR background) may be retrieved by a 
naive (and incorrect) assumption on the Lagrange multiplier fields in the Hamiltonian 
picture, or by dropping acceleration terms in the antisymmetric sector field equations 
that appear at higher-order, thereby revealing a constraint that does not properly 
exist in the theory. 

The relevant field configurations are those in which the antisymmetric components 
of the fundamental tensor may be considered as perturbations of dominant GR 
(symmetric sector) components on the initial data hypersurface Sq. Thus, an 
expansion in powers of the antisymmetric components (B a , ^[ ah l) about the symmetric 
components j( ab ) will be made, in which the components of the fundamental tensor 
take on the form: 



VS, (3 a ~ B a := S a bB b , Ct a ~ -J[ a b]B 



G(ab) ~ Sab, G[ab] ~ ~l[ab] '■= — S ac S ba ^ 



cd] 



(24) 



Terms of order n in (B a ,j^ ab ^) will be indicated by C(skew"), spatial indices will be 
raised and lowered using ^( ob ) and its inverse S ab , and '~' will indicate the dominant 
contribution in powers of (B a , <y afc l). 

The assumption is that all of the antisymmetric sector tensors (B a , W a , 7 ' ab ', j a b, 
W, b a: v a b , X bc ) are ©(skew 1 ), therefore leaving the dominant C(skew°) terms of the 
symmetric sector Lagrange multipliers identical to those of GR: 

L~0, c a ~0, a a ~^a b , u a b ~ k a b ~ ^ ac h bc , v£ 3 > [ 7 ] (ab) ~ 0, (25a) 
and those in the antisymmetric sector: 

W ~ IV a 3) [B] a , v% ~ V< 3) [fl] a - 15ZV^[B] C + j%, 
Kb - -h (cd) (vi 3) [ 7 ][ d a] + Vi 3 )[ 7 ] M - V( 3) [ 7 ] [ah] ) (256) 
+B a k b c - B b k a c - \5 c [a W b] + 5 c [a b b] . 

The trace of the last of these results in the constraints 

Xi w \l (ab) W b + V f ) [ 7 ] [afc] + kB a - k\B h w i 7 (Qb) TF fc + \ x S a - 0, (26) 

which are equivalent to (14) after dropping all contributions higher than ©(skew 1 ), 
and where the linear order contribution to S a in jl^ ) has been identified as (l/2) 1 S a ~ 

y(3)[ 7 ][a6] + kB a -k a b B b . 

Note that in dropping the contributions from the Lagrange multiplier b a , equation 
(14) has been required to play a very different role in the linearized system than in the 
nonperturbative case. Whereas in the exact treatment, the relation (14) is a condition 
determining the Lagrange multiplier b a , in this linearization it has become a constraint 
on the Cauchy data, and b a left at this stage as an undetermined Lagrange multiplier. 
As Xi are constraints, they must be preserved in time (to this order) for the evolution 
of the system to be consistent. It is straightforward to show that all contributions 
to xi involving b a cancel, resulting in an additional three constraints x% '■= Xi ~ 0- 
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Using the constraint ( |2q ) to determine W a , the linear contribution from (19a) clearly 
depends on b a , as does (20a), as it takes on the form 

W a « £f [W] a - 2N l Z [aX] - m 2 NB a , (27) 

where the linear order contribution to (216) is (R^ b is the O(skew ) GR Ricci tensor) 
l Z [aM ~ -R { a 3 b ] B b + k a b k b c B c - B a k b c k\ + k abl ^a c - 2 3ab a» 

+ 1 1 k\ a W b] 3k b [a b b] - (V< 8 > [B] b - 5 b a V& [Br)ab (28) 

-7 (bc M 3) b'] ac + V( 3 )[fc ab7 M] + fc hc vi 3) [ 7 ] [ca] . 

Considering (20) as second-order field equations for the configuration space 
variables, ( p6| ) determines b a as 

b a ~ \W a - \ lr E a ; (29) 

however, inserting this into (20a) yields the algebraic constraint (which is in fact 
equivalent to \2a ~ 0) 

W a ~ -% ~ £f [ X H] - 2N l Z [aX] ~ m 2 NB a , (30) 

instead of the expected evolution equation. (The field W a is easily determined in terms 
of B a using the linear contribution from (19a), and the time derivatives on k ab that 
occur in lr EL a are removed using (206).) The equations ( |30| ) then play a similar role as 
the time-space components of the field equations in the massive Kalb-Ramond theory: 
d a F 0aa + m 2 /i[° a l = 0, which do not involve acceleration terms and are therefore 
constraint equations. 

One could take this system seriously and determine whether the constraint algebra 
closes properly by computing X2 j presumably leading to a determination of b a instead 
of further constraints. This would recover a system analogous to the massive Kalb- 
Ramond field, and, consequently, the linear approximation of Moffat (1995a, b) and 
Clayton (1996a). A more systematic description has not been given since, as we 
shall see, this system is not an accurate representation of weak antisymmetric sector 
dynamics. It will be shown that no additional constraints properly exist, and the 
majority of choices of arbitrarily small antisymmetric sector Cauchy data result in 
singular evolution. Avoiding such configurations is possible by a particular choice of 
initial data, but in that case the system does not evolve in a manner consistent with 
the above linearization. 

Considering the exact second-order system, (19a) may be solved exactly for W a in 
terms of Cauchy data as 

W a = ~ ±=S ab (B b - £f[B] b + Ne c [ 7 }^ +N^C ±C ± ), (31) 

which, combined with (18d), may be used in order to determine the exact form of the 
acceleration terms appearing in (20 a) to be 

- ^j^=0- a \O bc S cd B d - (32) 

the time derivatives on k ab and j ab that appear are once again removed by making use 
of (206). Clearly to lowest order (9° b ~ j( ab \ and so its inverse is 0^ b ~ S a b- The 
operator Of 6 in (15a) on the other hand, is C(skew 2 ) 

Of - B a B b + -L 7 [ ac ] 7[cd]7 W = B a B b - 7^)7 • 7 + 7 a 7 b , (33a) 
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where (24) has been used in order to deduce that G a ~ — 7 a in (16), also leading to 
the form of the inverse from (156) 



1 B ■ B — 7- 7 i 

7(q6) - ^—^{Balb + laB b ) H ■ 7a7hJ , (336) 



7 • 7 

which is clearly 0(skew~ 2 ). The presence of the operator 0% b in ( |32] ) immediately 
shows how these acceleration terms disappear in the above linearized analysis, since 
it causes the acceleration terms in ([32]) to appear at third-order, leaving these field 
equations as constraints which do not truly exist in the theory. 

Furthermore ( |32| ) will result in very poorly behaved accelerations near 
configurations where the antisymmetric sector Cauchy data is small unless an ad hoc 
condition is imposed on the initial data, leaving a nontrivial wave equation for B a . 
To see this, (206) may be written to C(skew 3 ) as: 



l—L=0 2ab B b + % a ~ £fl 1 Z] a ~ 2N 1 Z [al _ ] - Nm 2 B a + 3 * Q , (34) 



where 3 ^> a are the remaining non-acceleration contributions to the field equations. 
The acceleration of B a is therefore given by 



B a ~ § WSO2 1 (-2 X2b + 3 * 6 ), (35) 



where \2a has been identified from the time evolution of ( pq ) . Due to the presence of 

O^ 1 " in (]35|), the right hand side is generally 0(skew _1 ), resulting in unboundedly 
large accelerations resulting from arbitrarily small antisymmetric sector Cauchy data 
on So. (Note that the sign of the right hand side of ( |35"| ) is uncorrelated with that of B a 
in general, since the sign of 0% b may be chosen by adjusting the relative magnitudes 
of B a and 7[ Qb l on E 

.) To avoid this, one could require that the initial data satisfy a 
condition that mimics the linearized result ( p6[ ) 

X2a ~ ty„, (36) 

effectively constraining the linear order contributions to (|35| ) to vanish up to terms of 
higher-order. However even in this case, (|35| ) becomes 

B Q ~ \NVS0 2 ia \^b + 3 *b), (37) 

which is a nontrivial ©(skew 1 ) evolution equation for B a that is not reproduced by 
the naive linearization. (Note that the second order spatial derivatives that combine 
with the time derivatives of B a to make (20a) a hyperbolic wave equation have also 
dropped out of Z[ a ±] at linear order in (28), presumably reappearing at third order 



to give (37) the appropriate hyperbolic form.) 

From the first-order point of view, the relation (186) determines the Lagrange 
multiplier b a in all cases when the operator is nondegenerate, including weak 
antisymmetric sector field configurations. (If 0% is degenerate then it may only be 
inverted on some subspace, determining some components of b a and leaving the rest 
undetermined. The analysis of these cases would be rather more complicated and will 
not be discussed further.) Since, as can be readily seen from (156), 0^ b is ©(skew -2 ) 
and therefore generally becomes singular for vanishingly small antisymmetric sector, 
the presence of b a in (19a) would drive B a to evolve arbitrarily quickly as one considers 
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vanishingly small antisymmetric sector Cauchy data. Explicitly, the last term in (19a) 
may be replaced using (186) leaving 

B Q « £f[B] a - iV(v£ 3) [7] [o6] +j [ab] C X b X ) + jN^O-^ c {OfW d + ~ c ). (38) 

The presence of in the last term indicates that unboundedly large velocities 

generally occur for infinitesimally small antisymmetric sector Cauchy data on So, and 
it is easy to see that the presence of b a in Z\ajJ will lead to similar behaviour for (20 a). 
Once again choosing initial data such that (g^) is satisfied removes all of these effects, 
as then (186) becomes (to 0(skew 3 )) 

b a ~ W* + k°2ab( 2 V - B b B c W c + 3 E b ), (39) 

effectively resulting in b a ~ C(skew 3 ). (Note that this alone does not guarantee 
that b a vanishes smoothly with the vanishing of the antisymmetric sector, but it 
is straightforward to choose data such that it does; for example, one could choose 
OfW b + S a ps Off b , where f b is an ©(skew) covector field. 

Although we have been focusing on the dynamics of (B a ,W a ), it is fairly 
straightforward to see that unless (|3£|) is imposed on Eo guaranteeing that b a ~ 
©(skew 1 ), the evolution of the symmetric sector will generally not resemble that of 
GR. The contributions from b a to the Hamiltonian constraint and the symmetric sector 
evolution equations imply that the symmetric sector cannot in general be considered 
as perturbative corrections to GR dynamics, since there are contributions from the 
antisymmetric sector that show up at 0(skew°). Even if one imposed ( |36| ) on £q, one 
would have to check that it is preserved in time in order for the system to remain well- 
behaved in evolution. This calculation would depend strongly on the chosen form of 3 ip 
and will not be pursued here. Nevertheless, it is clear that if one makes the choice (|36|), 
the resulting dynamics cannot be well approximated by the linearized system, simply 
because there are no additional constraints appearing. The evolution equation (19a) 
is not a constraint, and neither is (20a), showing that the full six degrees of freedom 
in the antisymmetric sector propagate even in the weak field regime, contradicting 
the linearized results. This behaviour in the antisymmetric sector indicates that the 
linearization given by Moffat (1995a, b) and Clayton (1996a) does not represent the 
weak-field evolution as determined from the full field equations, and therefore cannot 
be trusted: mNGT is not linearization stable about GR backgrounds. 

It is clear that these instabilities occur whenever the Lagrange multiplier b a 
becomes singular as the perturbations of the Cauchy data vanish. The general form 
of the solution of (14) specialized to spherically symmetric systems is 

^2W + ^' (40) 

and one finds that in regions of spacetime where B 1 is vanishingly small (for example, 
in an exact GR background or in the asymptotic region of an arbitrary asymptotically 
flat spacetime), perturbations generally cause b\ to become arbitrarily large, which 
in turn cause antisymmetric sector velocities to become large. Thus although the 
Wyman sector solution (Wyman 1950; Cornish 1994) (which is becoming the basis 
for much of the phenomenology of NGT (Moffat and Sokolov 1995a, b)) assumes that 
both B 1 and W\ vanish globally, if one considers perturbations of these fields on 
Y,q, the above behaviour reappears. Therefore, although the instabilities have been 
proven to exist for configurations very close to purely GR spacetimes, one expects 
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that any asymptotically flat spacetime will also suffer them. In particular, since there 
is no asymptotically flat spacetime in which (B 1 , W\) is asymptotically nonvanishing 
(Clayton 1996a), the unique asymptotically flat, static spherically symmetric solution 
with nontrivial antisymmetric sector (the Wyman solution) will also be unstable 
against perturbations. 

This is essentially the same effect as was found in (Isenberg and Nester 1977), 
where the effect of gravitational dynamics on the constraints of various derivative 
coupled vector fields was studied. It was found that constraints on the vector field 
may be lost when GR is considered as evolving concurrently with the vector field (as 
opposed to the vector field evolving on a GR background) . This manifested itself as an 
increase in the number of degrees of freedom in the vector field, and singular behaviour 
in the evolution equations when approaching asymptotically flat spacetimes. There 
one finds no evidence of this when considering the vector field dynamics on a fixed 
GR background, which is analogous to the linear approximation here. 

Kuchaf (1977) has noted that fields such as these that are derivatively coupled 
to GR (as NGT may be considered to be) may propagate off of the light cone as 
determined by the spacetime metric, and indeed, there are explicit examples of such 
behaviour (Buchdahl 1958, 1962; Cohen 1967; Velo and Zwanzinger 1969; Aragonc 
and Deser 1971). It is therefore not surprising that the field equations of UFT (and 
hence NGT) have been found to propagate information in a manner that is inconsistent 
with a simple causal structure as determined from a single spacetime metric; in fact 
Maurer-Tison (1959) has found that in UFT there are three such metrics that one 
must take into account. In UFT (or massless NGT), these metrics were consistent in 
the sense that there was locally a largest light cone (although the metric that defined 
this would change depending on the strength of the antisymmetric sector), and in 
the limit of vanishing antisymmetric sector all three degenerated to the single causal 
metric of GR (Maurer-Tison 1956). Although these metrics are not known for mNGT 
at this time, given the fact that the acceleration terms for B a have no simple GR 
limit, it is reasonable to suspect that one of these light cones has been altered so as 
to destroy the above mentioned consistency (perhaps resulting in a light cone that is 
degenerate in the GR limit). 

4. Conclusions 

It has been shown that arbitrarily small antisymmetric sector Cauchy data leads 
to singular evolution for the majority of possible choices of perturbatively small 
antisymmetric sector initial data. The results followed from an examination of the 
exact field equations of the massive nonsymmetric gravitational theory near purely 
general relativistic field configurations, and has been demonstrated by considering 
both the first and second order points of view, showing how the naive linearization 
fails to accurately describe the system, even in those cases where the choice of initial 
data does not lead to singular evolution. In doing so, it has been shown the constraints 
that guaranteed good fall-off for the linearized fields do not properly exist even for 
weak fields. Given this, it seems that the criterion for choosing the form of the mNGT 
action (Clayton 1996a) has not truly been fulfilled, and it remains unclear whether 
one has truly made an improvement over the original formulation of NGT. 

The failure of linearization stability has been noted in GR for closed spaces (Brill 
and Deser 1973). Here it has been found that not only can the linearized system not 
be trusted, but also that what appear to be benign perturbations of particular initial 
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data (i.e., as having very little effect on the evolution of the system as a whole) result 
in very different evolution, in which some of the antisymmetric sector fields are given 
arbitrarily large velocities. This is clearly interpretable as a Cauchy instability in the 
usual sense, since the evolution of generic configurations that are arbitrarily 'close' 
to a GR spacetime (or part thereof) does not smoothly depend on the initial data, 
and makes it difficult to physically interpret such spacetimes. Denying the physical 
importance of these configurations amounts to labelling as unphysical purely GR 
spacetimes and asymptotically flat spacetimes, and a Newtonian limit would instead 
have to be recovered in some (presumably stable) region of spacetime in which the 
antisymmetric sector is not small compared to the symmetric sector. Clearly if one 
does not allow asymptotically flat (or nearly so) spacetimes, nor regions of spacetime 
where the symmetric sector dominates to the point where one essentially recovers GR 
physics locally, then these instabilities are avoided. 

Despite the fact that these results indicate that the relevant (Wyman) solution 
would be unstable, there are some encouraging phenomenological results on galaxy 
dynamics (Moffat and Sokolov 1995a), as well as some optimism that the collapse of 
spherically symmetric matter would be nonsingular (Moffat and Sokolov 1995b). It is 
possible to further modify the dynamics of NGT in order to remove these instabilities 
and thereby recover this phenomenology by guaranteeing that three of the field 
equations appear as constraints rigorously, cither imposed via Lagrange multipliers 
in the action (Moffat 1996), or leading to generalizations of the model introduced 
by Damour et al (1993) that possess 'gauge invariant kinetic terms'. Alternatively, 
given that the antisymmetric sector should appear as a Kalb-Ramond field (Kalb and 
Ramond 1974), it is reasonable to expect that there should be some tie between NGT 
and string theory (Moffat 1995d). It may then be the case that the fundamentally 
nonlocal nature of strings (which does not have a Cauchy initial value formulation) 
would cause the NGT field equations to be an inadequate description of the system 
in some situations. 

Acknowledgments 

I would like to thank my supervisor J. W. Moffat and the University of Toronto for 
supporting this research, J. Legare for his careful reading this manuscript, and J. W. 
Moffat and L. Demopoulos for comments and discussions related to this work. 

References 

Aragone C and Deser S 1971 II Nuovo Cimcnto 3A 709-20 
Arnowitt R, Deser S and Misner C W 1959 Phys. Rev. 116 1322-30 
Bonnor W B 1954 Proc. R. Soc. London A226 366-77 
Brill D and Deser S 1973 Comm. Math. Phys. 32 291-304 
Buchdahl H A 1958 II Nuovo Cimcnto 10 96-103 

—1962 II Nuovo Cimcnto 25 486-96 
Callaway J 1953 Phys. Rev. 92 1567-70 

Clayton M A 1995 Hamiltonian formulation of nonsymmetric gravitational theories University of 

Toronto Preprint UTPT-95-20 
—1996a J. Math. Phys. 37 395-420 
1996b Massive Nonsymmetric Gravitational Theory: A Hamiltonian Approach Ph.D. thesis, 

Department of Physics, University of Toronto 
Cohen H A 1967 II Nuovo Cimento 52A 1242-52 
Cornish N J 1994 Mod. Phys. Lett. A 9 3629-40 

Damour T, Deser S and McCarthy J 1992 Phys. Rev. D 45 R3289-91 



14 



—1993 Phys. Rev. D 47 1541-56 
Einstein A 1945 Ann. Math. 46 578-84 
Einstein A and Straus E G 1946 Ann. Math. 47 731-41 

Hawking S W and Ellis G F R 1973 The Large Scale Structure of Space- Time (Cambridge: Cambridge 

Monographs on Mathematical Physics, Cambridge University Press) 
Infeld L 1950 Acta. Phys. Pol. 10 284-93 
Isenberg J A and Nestor J M 1977 Ann. Phys. 107 56-81 

— 1980 General Relativity and Gravitation: One Hundred Years After the Birth of Albert Einstein 

vol I cd Held A (New York: Plenum Press) p 23-97 
Kalb M and Ramond P 1974 Phys. Rev. D 9 2273-84 
Kelly P F 1991 Class. Quantum Grav. 8 1217-29 

1992 Class. Quantum Grav. 9 1423, Erratum 

Kuchaf K 1977 J. Math. Phys. 18 1589-97 

Kunstatter G, Leivo H P and Savaria P 1984 Class. Quantum Grav. 1 7—13 
Legare J and Moffat J W 1995 Gen. Rel. Grav. 27 761-75 
Maurer-Tison F 1956 C. R. Acad. Sc. 242 1127-9 

1959 Ann. Scient. Ec. Norm. Sup. 76 185-269 

Moffat J W 1979 Phys. Rev. D 19 3554-58 
—1984 Found. Phys. 14 1217-52 

— 1991 Gravitation: A Banff Summer Institute ed Mann R and Wesson P (New Jersey: World 

Scientific) p 523-97 
—1995a Phys. Lett. B 355 447-52 

1995b J. Math. Phys. 36 3722-32 

1995c J. Math. Phys. 36 5897-915 

— 1995d Nonsymmetric gravitational theory as a string theory University of Toronto Preprint 
UTPT-95-26 

— 1996 Dynamical constraints in the nonsymmetric gravitational theory University of Toronto 
Preprint UTPT-96-05 

Moffat J W and Sokolov I Yu 1995a Galaxy dynamics in the nonsymmetric gravitational theory To 
appear in Phys. Lett. B 

— 1995b On gravitational collapse in the nonsymmetric gravitational theory University of Toronto 
Preprint UTPT-95-21 
Teitelboim C 1973 Ann. Phys. 79 542-57 

— 1980 General Relativity and Gravitation: One Hundred Years After the Birth of Albert Einstein 
vol I cd Held A (New York: Plenum Press) p 195-225 
Vclo G and Zwanzinger D 1969 Phys. Rev. 188 2218-22 
Wyman M 1950 Can. J. Math. 2 427-39 



